To investigate the nature Ml-78 further, we present a de- 
tailed spectroscopic study of Ml-78 in the optical and near- 
infrared. We obtained long-slit, intermediate-resolution, op- 
tical spectroscopy with the ISIS spectrograph mounted on 
the William Herschel Telescope (WHT) at Roque de los 
Muchachos Observatory (La Palma, Spain). As a complement, 
we obtained long-slit, intermediate-resolution, near-infrared 
spectra using LIRIS, the near-infrared imager/spectrographer 
also installed at the WHT. Ml-78 is a high-density nebula with 
substantial physical differences between its two main morpho- 
logical zones: a bright arc to the SW and a blob of emission in 
the NE. Specifically, the blob in the NE has a higher electron 
QQ I temperature (13400 K) and visual extinction (about 9 mag) than 
' the SW arc. The most important result, however, is the confir- 
mation of a nitrogen enrichment in Ml-78. This enrichment is 
stronger at the location of the NE blob and is correlated with a 
defficiency in the O abundance and a (dubious) He enrichment. 
Such an abundance pattern is typical of ejecta nebulae around 
evolved massive stars such as Wolf-Rayet and Luminous Blue 
Variable stars. The spatial variations in the physical condi- 
tions and chemical abundances and the presence of more than 
one possible ionizing source indicates, however, that Ml-78 
, is better described as a combination of a compact H ii region 
+ ejecta. This is confirmed by the Hei2.112yum/Br7 line ra- 
Q ! tio, which indicates a hot (1^// > 40000 K) O star in the SW 
' arc. Finally, we detect H2 emission that extends over a large 
, (~ 30") area around the ionized nebula. Analysis of the near- 
infrared H2 lines indicates that the excitation mechanism is UV 
fluorescence. 
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Abstract. There is considerable controversy surrounding the nature of Ml-78, a compact nebula located beyond the Perseus 
arm. It was first classified as a planetary nebula and is nowadays generally considered to be a compact H ii region. 



1. Introduction 

Ml-78 (ano ther designa t ion is IRAS 21 190+5140) was first 
observed by 'Minkowski] ( 1946h in a survey of emission-line 
objects and is nowadays generally classified as an Hii re- 
gion. Nevertheless, it is still listed in the Simbad Astronomical 
Databas^ as a "possible planetary nebula" (PN), which re- 
flects the controvery that has existed aro und the nature of this 
source. The survey of Latter et alJ ( 1995 ) catalogues this object 
as a middle-butterfly PN, where a butterfly or bipolar nebula 
is a PN that is bi-lobed and has an equatorial waist. Phillip^ 
(|2004 places this source at a distance of 1.6 kpc using the 
correlation between 5 GHz brightness temperatures and intrin- 
sic radio luminosities established for PNe and based on stel- 
lar evolutionary mo dels. However, the distance determined by 
Puche et aU jl988h based on the Hi absorption spectrum is 
in between 5 and 8 kpc, beyond both the Local arm and the 

Perseus arm. This is in agreement with the kinematic distance 

I ■ 1 

of 8.9 kpc obtained from its radial velocity (e.g.'Peeters et al., 
20021). Such a distance implies a h i gh lu minosity of about 
2 X 10^ L/Lq (see again IPeeters et al.ll2002h . which makes the 
classification of Ml-7 8 as a PN untenable, as pointed out al- 
ready by Z iilstra et al.l ( 1990) and Gussie ( 1995). A classifica- 
tion as an ultracompac t H ii region is favored instead by these 
authors. [Gussie ( 1995 ) even favors a rather unique classifica- 
tion for this region, that of an ultracompact H ii region with a 
post-main sequence central star, probably a Wolf-Rayet (WR) 
star, based on the high expansion velocity of the nebula (around 
25 kms~' ) and a possible overabundance of nitrogen (see also 



Perinottd.il99lh . 

Several near-infrared (NIR) spectra of Ml-78 have been 
obtained showing emission of atornic hydrogen and helium 
and mo lecular hydrogen (e.g. fHora et aliri999HLumsden et all 
2001allbl) . This source has also been imaged in the NIR (e.g. 



Send offprint requests to: N.L. Martin-Hernandez, leticia@iac.es 
' http://simbad.u-strasbg.fr/simbad/ 



Latter et aU 1995 ) mid- infrared (e.g. Peetersl 2002 ) and radio 
(e.g.l Ziilstra et al. . 1 19901) . showing a similar morphology at all 
wavelengths, i.e. an overall structure emcompassing a promi- 
nent arc of which the tail at both sides converges towards a blob 
of emission to the NE (see Fig. [1). The ISO spectrum of this 
object is rathe r unique and shows silicate in emission at about 
10 and 20 //m ( iPeeters et al. [ I2OO2I) . In addition, it shows strong 
emission features generally attibuted to emission by polyclyclic 
aroniatic hydrocarbon (PAH) molecules, typical of H 11 regions 
(cf. |Peetersetal.[i2002l) . 

We obtained long-slit, intermediate-resolution, optical 
spectroscopy with the ISIS spectrograph mounted on the 
4.2 m William Herschel Telescope (WHT) at Roque de los 
Muchachos Observatory (La Palma, Spain). As a comple- 
ment, we obtained long-slit, intermediate-resolution, NIR spec- 
tra of Ml-78 using LIRIS, the NIR imager/spectrographer also 
mounted on the WHT. The aim of this study is to investigate 
the controversial nature of Ml-78. 

This paper is structured as follows. Section|2]describes the 
observations and data reduction. Section |3] and Sect. |4]present 
the results of the optical and NIR data, respectively. Section |5] 
discusses the nature of Ml-78 and Sect.|6]summarizes the main 
conclusions. 



2. Observations and data reduction 

2.1. Optical observations 

Intermediate-resolution spectroscopy was obtained on 2007 
July 18 with the ISIS spectrograph mounted on the 
4.2 m William Herschel Telescope (WHT) at Roque de los 
Muchachos Observatory (La Palma, Spain). Two different 
CCDs were used at the blue and red arms of the spectrograph: 
an EEV CCD with a configuration 4100 x 2048 pixels with a 
pixel size of 13.5 i^m in the blue arm and a REDPLUS CCD 
with 4096 X 2048 pixels with a pixel size of 15 yum in the red 
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arm. The dichroic used to separate the blue and red beams was 
set at 5300 A. The slit was 3.'7 long and I'.'OS wide. Two grat- 
ings were used, the R1200B in the blue arm and the R316R in 
the red arm. These gratings give reciprocal dispersions of 17 
and 62 A mm and effective spectral resolutions of 0.86 and 
3.56 A for the blue and red arms, respectively. The blue spectra 
cover from AA4225 to 5075 A and the red ones from /1/15335 to 
8090 A. The spatial scale is 0'.'20 pixel ' and a.'22 pixel ' in 
the blue and red arms, respectively. The average seeing during 
the observations was ~0'.'8. 

We observed a single slit position centered on the brightest 
nebular spot (see Fig.[Tll and with a position angle RA.=50°. 
We took three consecutive 600 s exposures in both arms in 
order to get a reasonably good spectrum of this faint object. 
The data were wavelength calibrated with a CuNe+CuAr lamp. 
The correction for atmospheric extinction was performed us- 
ing the average curve for continuous atmospheric extinction at 
Roque de los Muchachos Observatory. The absolute flux cal- 
ibration was achieved by observations of the standard stars 
BD-H33 2642, BD-H28 4211 and BD-i-25 4655. AH the CCD 
frames were reduced using the standard IRAfU TWODSPEC 
reduction package to perform bias correction, flat-fielding, 
cosmic-ray rejection, wavelength and flux calibration, and sky 
subtraction. 

2.2. NIR observations 

NIR spectra in the 2.0-2.4 fim wavelength range (/T-band) were 
obtained on the night of 2006 September 9 at the 4.2 m WHT 
using LIR IS, a long-slit, interrnediate - resolution, infrared spec - 



trograph (lAcosta Pulido et al.L 120031: iManchado et al.L |2004|) 



The spectral region blueward of about 1.98/zm is affected 
by substantially degraded atmospheric transmission. LIRIS is 
equipped with a Rockwell Hawaii 1024 x 1024 HgCdTe array 
detector The spatial scale is 0'.'25 pixel ', and the slit width 
used during the observations was 0'.'75, allowing a spectral 
resolution R - A/AA ~ 1000. The slit was orientated along 
RA.= 50° (slit #1, which corresponds to the single-slit ob- 
served in the optical) and 140° (slit #2) centered on the bright- 
est nebular spot (see Fig.[T]i. Weather conditions were relatively 
good, although with sparse cirrus. The seeing during our obser- 
vations varied between 0'.'8 and 1'.'3 measured from the FWHM 
of the standard star spectra. 

Observations were performed following an ABBA 
telescope-nodding pattern, placing the source in two positions 
along the slit, separated by 45". Individual frames were taken 
with integration times of 300 s and total on-source integration 
times of 30 minutes for slit #1 and 20 minutes for slit #2. In 
order to obtain the telluric correction and the flux calibration, 
the nearby HIP 106393 AO V star was observed with the same 
configuration. The data were reduced following standard 
procedures for NIR spectroscopy, using IRAF and LIRISDR, 
the LIRIS Data Reduction package. After the flat-field cor- 




^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
the Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation (http://iraf.noao.edu/ 1. 
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Fig. 1. NIR broad-band (Ks filter) acquisition image of Ml- 
78. The positions of the NIR slits (in red slit #1 with a P.A. 
of 50°, and in cyan slit #2 with a P.A. of 140°) are marked. 
The (0",0") position coiTesponds to R.A.=2l''20'"44.91 ' and 
Dec.=-H51''53'"27.3' (J2000.0). The Hj regions analysed in 
Sect. |4.7| are indicated. At a distance of 8.9 kpc, 1" corresponds 
to 0.043 pc. 



rection, consecutive pairs of AB two-dimensional spectra 
were subtracted to remove the sky background. The resulting 
frames were then co-added to provide the final spectrum. Sky 
lines present in the science data were used to determine the 
wavelength calibration. The final uncertainty in the wavelength 
calibration is roughly lA. The resulting wavelength-calibrated 
spectra were divided by a composite spectrum to remove 
the telluric contamination. This composite spectrum was 
generated from the observed spectra of HIP 106393, divided 
by a Vega model convolved with the actual spectral resolution 
(~ 23A) using the routine xtellcor^eneral within the Spextool 
packag^ an IDL-based spectral reduction tool written by 
Mike Cushing and Bill Vacca. This routine is also used for 
the flux calibration, which is done by normalizing with the B 
and V magnitudes of the standard star provided by the Simbad 
Astronomical Database (B=5.734 and V=5.766). We estimate 
an uncertainty in the flux calibration around 40% based on the 
comparison of the various spectra obtained for the standard 
star. 



3. Results from optical observations 

3.1. Features of the optical spectra 

The two-dimensional optical spectrum shows two peaks in the 
surface brightness distribution of the emission lines along the 
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Fig. 2. Red range of the optical spectra extracted from a slit 
position coincident with slit #1 of the NIR observations (see 
Fig. [1). The upper panel corresponds to the emission of the NE 
blob; the inset show an enlargment of the spectral area around 
the [Nil] 5755 A and Hei 5876 A hnes. The middle and the 
lower panels show the red spectrum of the SW arc at different 
line-flux levels. The fluxes are not corrected for reddening. The 
identification of most emission lines and absorption features 
are indicated. The broad absorption features longward 6800 A 
are due to the atmosphere. 

spatial direction. These peaks coincide with the two main mor- 
phological zones of the nebula: the brightest arc to the SW and 
the blob of emission to the NE. In addition, each peak con- 
tains a featureless, very weak and unresolved stellar continuum 
(note that the nebular contribution to the continuum in the op- 
tical is generally very small). Only the continuum associated 
with the S W arc - brightest one - is detected in the spectrum of 
the blue arm. We extracted two one-dimensional spectra isolat- 
ing the emission from each peak, the spatial extension of each 
zone being 4'.'4. The spectra of the red range of both zones are 
shown in Fig.|2] A detailed inspection of the emission line pro- 
files along the slit does not show any evidence of stellar features 
or broad components in the recombination-line profiles due to 
the contribution of photospheric emission of massive stars. 

Line intensities were measured integrating all the flux in 
the line between two given limits and over a local continuum 
estimated by eye and using the SPLOT routine of the IRAF 
package. The line intensities of the blue arm spectra were nor- 
malized to that of HyS, while those of the red arm were nor- 
maUzed to Ha. The reddening coefficient, c(Hj3), was deter- 
mined by fitting the observed /(Hy)//(Hy6) ratio of the spec- 
tra of the blue arm to the theoretical Hne rati o obtained from 
the calculations of IStorev & Hummeii (1 1995b . assuming case 



Once the reddening was derived, the emission line ratios of the 
red spectra were re-scaled to the 1(11/3). The spectrum of the 
NE blob shows a somewhat larger reddening coefficient, indi- 
cating a larger amount of internal dust in this zone of Ml-78 
that produces an extinction Ay of about 9 mag. Our determi- 
nations of cjHB) are consistent with the valu e of 4. 1 found by 



Aller & KevesI (Il987h . iLumsden et al.l (1200 lb ) derived a simi- 
lar extinction in the NIR of Ay ~ 6.5 mag (assuming Ry = 3; 



MathisL 11990) estimated by comparing the higher Brackett se- 



ries lines with Bry and assuming an exponential extinction law 

In Table [T] we include, for the two extracted areas, the ob- 
served and laboratory wavelength of the identified lines, their 
reddening-corrected intensity ratios with respect to H/3, the cor- 
responding values of the reddening function, f(A), for each 
emision line, the reddening coefficient c(Hy6), the observed H/S 
line flux and the physical conditions of the gas - electron den- 
sity and temperature - determined from different emission-line 
ratios. We performed a single Gaussian fit of the profiles of 
the brightest lines of the blue range of the spectra (H/3 and 
[O m] 4959 and 5007 A) for determining the radial velocity of 
the nebula with respect to the local standard of rest (LSR), 
finding a value of -74+1 and -70±2 km s"' for the SW arc 
and the NE blob, respectively. These values a re in good agree- 
ment w ith the Vlsr - -76+2 km s ' found bv lGussie & Taylor 
(1 1989b also from optical spectroscopic observations and with 
m easurements o f hydro gen radio recombination line s obtained 
bv lTerzianelal] ( Il974 and IChurch well et all (Il976b . 

The red spectrum of the NE blob shows several uncommon 
permitted faint emission lines that were identified as O i 7254 
A (a blend of three unresolved lines at 7254.14, 7254.45 and 
7254.53 A), 1 7002 A (a blend of 7001 .92 and 7002.23 A) and 
Oi 6046 A (a blend of 6046.23, 6046.44 and 6046.49 A), be- 
longing to multiplets no. 20, 21 and 22, respectively; Ni 7442 
and 7468 A of multiplet no. 3; and Siii 6347 and 6371 A of 
multiplet no. 3. As demonstrated by iGrandil (119751 Il976b . all 
these lines come from upper levels of their respective ions that 
can be excited from the ground state by the absorption of stellar 
photons longward of 912 A. Also, it is important to note that 
all these lines are precisely the bright est permitted ones of O i. 



N I and Si ii in the Orion Nebula fsee lEsteban et al.L l2004^ . al- 
though their line ratios with respect to H/3 are larger in Ml-78 
than in the Orion Nebula. These permitted lines are produced 
in the vicinity of the transition zone of the nebula, where the 
fraction of neutral mater ial becomes im portant and are rather 
uncommnon in PN (e.g. iLiu et al.l 12004) . where they should 



be presumably produced in dense clumps or g lobules of neu - 
tral material embedded inside the ionized gas (" Grandi . 1976 ). 
Other interesting lines in the spectra of Ml-78 are those of 
[Fell] 7155 A (multiplet 14F, only detected in the NE blob) 
and [Ni ii] 7378 A (multiplet 2F, detected in both spectra). Both 
are also the brightest forbidden lines of those ions in the Orion 
Nebula in the s pectra l range covere d by our red arm spectrurn 



(lEsteban et al. 



2004 . ILucvI d 1995b and iBautista et all d 1996b 



B, and electron temperature and density of 10"^ K and - 10 ^ 
cm"^. We considered the reddening law of iWhitfordI (1958). 



demonstrate that, in the Orion Nebula, these lines are pro- 
duced by the absorption of stellar non-ionizing UV photons. 
Interestingly, these lines, especially those of [Niii], are very 
strong in the circumstellar ejecta around the Luminous Blue 
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Table 1. Dereddened line intensity ratios with respect to 1(11/3) =100, reddening coeeficient, observed H/3 flux and physical 
conditions obtained from the optical spectra. 



SWarc 



NE blob 





Identiiication 






1(A)/ 1 (tip) 




l(A)/l(tip) 


4340.47 


Hy 


0.135 


4339.00 


47.0 ± 2.5 


4339.04 


47.0 ± 3.9 


4471.48 


He I 


0.10 


4470.09 


4.0 ± 0.8 






4861.33 


UJS 


0.00 


4859.73 


100 ± 2 


4859.77 


100 ±2 


4921.93 


He I 


-0.01 


4920.42 


9.2: 






4958.91 


[Om] 


-0.02 


4957.30 


134.6 ± 3.8 


4957.38 


93.6 ± 2.7 


5006.84 


[Om] 


-0.03 


5005.17 


372.1 ± 8.3 


5005.26 


258.3 ±5.8 


5015.68 


He I 


-0.03 


5013.99 


2.1 ±0.4 


5014.15 


2.0 ± 0.4 


5517.71 


[CI III] 


-0.17 






5516.19 


0.8 ± 0.2 


5537.88 


[Clm] 


-0.18 






5536.28 


1.0 ±0.2 


5754.64 


[Nil] 


-0.21 


5752.77 


1.3 ±0.2 


5752.69 


2.5 ±0.1 


5875.67 


He I 


-0.23 


5873.75 


11.4 ±0.4 


5873.72 


19.8 ± 0.6 


6046.44 


Oi 


-0.26 






6044.87 


2.5: 


6300.30 


[Oi] 


-0.30 


6298.55 


0.9 ±0.1 


6298.50 


2.4 ±0.1 


6312.10 


[Sra] 


-0.30 


6310.38 


1.8 ±0.1 


6310.27 


2.4 ±0.1 


6363.78 


[Oi] 


-0.31 


6362.32 


0.2: 


6361.96 


0.7 ±0.1 


6571.36 


Sin 


-0.31 






6370.07 


2.9 0.7 ± 0.4 


6548.03 


[Nil] 


-0.34 


6546.10 


18.5 ±0.5 


6546.02 


26.4 ± 0.8 


6562.82 


Ha 


-0.34 


6560.85 


286.0 ± 5.9 


6560.74 


286.0 ± 5.9 


6583.41 


[Nil] 


-0.34 


6581.41 


59.4 ± 1.3 


6581.32 


59.4 ± 1.3 


6678.15 


He I 


-0.35 


6676.08 


3.8 ±0.1 


6676.00 


3.8 ±0.1 


6716.47 


[Sii] 


-0.36 


6714.33 


2.1 ± 0.1 


6714.05 


2.1 ±0.1 


6730.85 


[Sii] 


-0.36 


6728.76 


3.8 ±0.1 


6728.46 


3.8 ±0.1 


7001.92 


Oi 


-0.39 






6999.58 


0.2: 


7065.28 


He I 


-0.40 


7062.74 


8.8 ±0.3 


7062.70 


8.8 ± 0.3 


7135.78 


[Arm] 


-0.41 


7133.23 


17.4 ±0.5 


7133.11 


17.4 ± 0.5 


71 1 A 










71 S9 9S 


n 1 • 




n I 


f) d.9 






7251.66 


2.3 ± 0.8 


/ ZoL.JJ 


rle I 


-U.4j 


7279.25 


1.0 ±0.2 


7278.96 


0.43 ± 0.04 


7^18 '^Q 

/ J i O. J^' 


IIJ 


-0 4^ 


7317.43 


9.8 ± 0.4 


7317.18 


8.7 ± 0.3 


7329.66 


[Onl 


-0.43 


7327.83 


9.1 ±0.3 


7327.59 


7.5 ± 0.2 


7377.83 


[Niii] 


-0.44 


7375.68 


0.2: 


7374.98 


0.20 ± 0.02 


7442.30 


Ni 


-0.44 






7438.75 


0.09 ± 0.02 


7468.31 


Ni 


-0.45 






7465.74 


0.15 ± 0.03 


7751.10 


[Arm] 


-0.48 


7748.92 


4.5 ±0.1 


7748.79 


1.56 ±0.06 


c(m 






3.09 


±0.02 


4.10 


±0.04 








160 + 3 


73 ±2 


«,([Sii])W 






6700 


± 1000 


6400 


± 1200 


ne([Cliii])« 










54Q0+8800 


7'e([Nll])<''> 






10900 + 600 


13400 ±400 



Blend of Oi lines at 7254.14, 7254.45 and 7254.53 A. Observed (uncorrected) flux in units of 10" 16 erg cm"^ s"'. 



In cm-\ In K. 



Variable (LBV) star P Cygni d johnson et"an. 1 1992 1) and an- 
other two nebulositie s associated with LBV stars in the LMC 
dStahl & WolA 19861) . Stellar fluorescence also seems to be the 
most likely excitation mechanism of [Ni n] lines in those nebu- 
lae (Lucy , 1995). 

As it can be seen in Fig. |2] both one-dimensional spectra 
show clear absorption features. Although the most conspicuous 
ones are atmospheric bands, the narrower features are of inter- 
steUar origin. We identified the Nai doublet at 5890 and 5896 
A, the K I doublet at 7665 and 7699 A, and several diffuse inter- 
stellar bands (DIBs) at 5780, 5797, 6177, 6284, 6379 and 6614 



A. We determined the LSR velocities of the centroids of the 
single Gaussian fits of the interstellar absorption features com- 
paring the observed and rest frame wavelengths (for the DIBs , 



we used the precise wavelengths of lGalazutdinov et alil2000h . 
The vlsr values we found are in the range -68 to -32 kms~\ 
consistent with the velocities of the H i absorption features ob- 
served bv I Puche et al. ( 1988 ) and identified as caused by neu- 
tral gas in and beyond the Perseus spiral arm. The relatively 
high negative radial velocity of the nebula and the presence of 
these interstellar absorption systems in the line of sight of the 
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Table 2. Ionic and total abundances derived from optical lines. 



SWarc 



NE blob 



12+log(0+/H+) 


7.81 ± 


0.17 


7.31 


+ 


0.10 


12+log(0++/H+) 


8.00 ± 


0.09 


7.59 


+ 


0.05 


12+log(0/H) 


8.24 ± 


0.10 


7.77 


+ 


0.06 


12+log(N+/H+) 


7.00 ± 


0.07 


6.93 


+ 


0.04 


log(N/0) 


-0.81 H 


t0.21 


-0.38 d 


:0.12 


12+log(S+/H+) 


5.43 ± 


0.10 


5.52 


+ 


0.07 


12+log(S++/H+) 


6.47 ± 


0.12 


6.27 


+ 


0.06 


12+log(S/H) 


6.53 ± 


0.12 


6.39 


± 


0.06 


log(S/0) 


-1.69 H 


t0.14 


-1.38 d 


:0.09 


12+log(Ar++/H+) 


6.11 ± 


0.07 


5.70 


+ 


0.05 


12+log(Cl++/H+) 






4.78 


± 


0.13 


He+/H+ (4471) 


0.078 ± 


0.019 








He+/H+ (5876) 


0.077 ± 


0.005 


0.119 


+ 


0.006 


He+/H+ (6678) 


0.097 ± 


0.006 


0.083 


+ 


0.004 


He+/H+ (7065) 


0.175 ± 


0.019 


0.079 


+ 


0.010 


<He+/H+><'"" 


0.084 ± 


0.009 


0.101 


+ 


0.018 


He/H 


0.091 ± 


0.009 


0.123 


+ 


0.018 



Mean value excluding He+/H+ (7065 A). 

object clearly indicate that Ml-78 is a rather distant Galactic 
object. 

3.2. Physical conditions and chemical abundances 

The physical conditions of the nebula, electron density and 
temperature were derived with the IRAF task TEMDEN of the 
pa ckage NEBULAR based in the FIVEL pr ogram developed 



bviP e Rober tis et al.l (Il987h and improved by Shaw & Dufoui 
( 11995) . The electron density, ris, was calculated from the [S ii] 
6717/6731 line ratio in the two one-dimensional spectra as 
well as the [Clm] 5518/5538 line ratio in the spectrum of 
the NE blob. The densities found are very high and simi- 
lar in both zones, 6700 and 6400 cm"^ for the arc and the 
blob, respectively (see Table [1]), and are roughly consistent 
but slightly higher t han previous determin ations based on opti- 



cal (~ 2500 cm ^ lAller&Keve"sill987h or far-infrare d spec 



troscopy (~ 3000 cm'^.'Ma rtm-Hernandez et al.Ll2002h . These 
large values of We are common in compact H ii regions and are 
also on the order of those observed in the Orion Nebula (e.g. 
Estebanet all ll998V 



Although the excitation degree of Ml-78 is rather high and 
the optical nebular [Oiii] lines are strong, the auroral [Oiii] 
4363 A line cannot be detected due to the large dust extinc- 
tion. Fortunately, the temperature sensitive [Nn] 5755 A line 
was measured with a rather good signal-to-noise ratio in both 
zones and therefore we could obtain a direct determination of 
the electron temperature, T^., in Ml-78. The values of Tg are 
included in Table [1] It can be seen that the temperature is sub- 
stantially different in the two zones; the NE blob shows T^. = 
13400 K, which is atypically high for a Galactic Hii region. 
A previous det ermination of the temperature was obtained by 
Aller&Keyesl (fl987). who found Te = 11 300 K - probably de- 
termined from the [O ii] line ratio - in agreement with the val- 



ues we obtain for the S W arc. Since the intensity of the [N ii] 
5755 A line could suffer from some contribution due to recom- 
bination, w e estimated this contribution using the expression 
obtained by Liu et al. ( 2000h and find that it is negligible, ap- 
proximately 0.4% and 0.2% for the arc and the blob zones, re- 
spectively. 

Ionic abundances of N+, 0+, 0++, S+, S++, Ar++ and C1++ 
were derived from collisionally excited lines by making use of 
the IRAF task IONIC of the package NEBULAR. In the ab- 
sence of temperature determinations for further ions, we con- 
sidered a one-zone scheme, assuming that the temperature ob- 
tained from the [N ii] line ratio is valid for the rest of the ions in 
the nebula. We measured several He i lines in our spectra. These 
lines arise mainly from recombination but they can be affected 
by collisonal excitation and self-absorption effects. We derived 
the He^/ H^ ratio using the effect ive recombination coefficients 
given by Pequignot et alJ j 1991 ) and applying the collisional 
corrections proposed by Kin gdon & Ferlandl (Il995h . The final 
adopted value for the He^ abundance is the mean of the values 
obtained for the brightest individual He i lines but excluding 
He I 7065 A, which suffers from the largest collisional effects. 
All the ionic abundances are included in Table |2] 

Since we do not detect He n lines, we expect that the contri- 
bution of O^^ should be negligible, and the total O abundance 
was therefore determined by just adding the ionic abundances 
of O^ and O^^. As can be seen in Table |2] the O/H ratio ob- 
tained for the two zones of the nebula is very different, the SW 
arc shows a value consistent with that expected by the ra dial 
Galactic abundance gradient of e.g. lDeharveng et al. I (l2000l) for 
a galactocentric distance of 12.7 kpc. However, the spectrum of 
the NE blo b shows an O/H ratio a factor of 3 lower than that 
of the arc. lAUer & Key e si (Il987n also find a relatively low O 
abundance for the object (12+log(0/H) = 8.07). 

The observed N^/0^ ratio - that is usually assumed to be 
equal to the N/O ratio due to the similar ionization potential 
of both ions - is high in both areas of the nebula, especially 
in the blob, where log(NVO^) - -0.38. The expected value of 
log( N/0) at the distanc e of the nebula should be about -1.2 



(e.g. Shaver et al.Lll983l) . indicating that the nebula is nitrogen- 
rich by factors as large as 6 in the NE blob. The presence of 
nitrogen enrich ment in the object ha s bee n indicated by se veral 
authors, such as lZijlstra et alJ (Il990l) and lPerinottol (Il99ll) . 

We also estimated the total abundances of He and S by 
making use of an appropriate ionization correction factor, icf, 
scheme for each element. These factors account for the contri- 
bution of the abundances of unseen ionization stages of each 
ch emical element. In the case of He, we used the /c/ proposed 
by Peimb ert et al.l (119921) . based on the similarity of the ion- 
ization potentials of He" and S^. This /c/ seems to work well 
for relatively high excitation objects like Ml-78, where only 
a small fraction of He" is expected. In the case of S, we have 
to consider some contribut ion for unseen S ^^ and we assume 
the correction proposed bvlStasinska (Il978h . The e quations of 
both icfs can be found in lGarcia-Roias et alj ( 2004 V As shown 
in Table |2] the S/H ratios of the two one-dimensional spectra 
differ by only 40% (0.14 dex) and can be considered rather 
similar within the errors. The fact that the total abundance of S 
does not show the behavior of the O/H ratio seems to support 
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that the difference of O abundance in both zones of the nebula 
could be real and not an artefact produced by, for example, a 
localized abnormally high T,.. On the other hand, the He/H ra- 
tio seems to be somewhat higher (34%) in the NE blob than in 
the SW arc. However, this difference is produced entirely by 
the different He^/H^ ratio obtained from the He i 5876 A line 
in both spectra and is not correlated with similar differences in 
the abundance obtained from the other He i lines. Therefore, we 
consider the different He/H to be dubious. 

4. Results from NIR observations 

4.1. Features of the NIR spectra 

Figure [3] shows the integrated nebular spectra of Ml-78 along 
slits #1 and #2. The most obvious features in the spectra are 
the Bry line, a large set of H i Pfund lines, several He i emis- 
sion lines and a relatively large number of molecular hydrogen 
transitions. Three forbidden lines of [Fern] are also evident. 
Table|3]lists the observed lines, their identifications and fluxes. 
The quoted uncertainties are those coming from the line fitting 
procedure. Peak positions, equivalent widths and fluxes of the 
lines were measure by fitting a Gaussian profile. A line was 
defined as being detected if its peak intensity exceeds the rms 
noise level of the local continuum by at least a factor of 3. 

The validity of our flux calibration can be checked by 
comparing the value we obtain for the Bry line flux with 



Table 3. NIR line fluxes extracted across slit #1 (between -7'.'5 
and 0") and slit #2 (between -2" and 2"). 



19991) 



previous spectroscopic observations. iHora et al 
sured a flux of 1110 x 10"^ W m"^ with a l"x6" slit ori- 
entated north-south and cen tered on the br ightest infrared por- 
tion of the source. Similarly. iLumsden et al. (2001b) measured 
1430 X 10"''^ W m"- integrating over a region l'.'2xl2'.'6. These 
values are comparable to the flux we observe for Bry and list 
in Tabled 

4.2. The spatial variation of NIR emission lines 

Figure |4] illustrates the spatial variation of Bry, the strongest 
He I lines (at 2.058 and 2.112 yum) and the H2 1-0 S(l) tran- 
sition at 2.121 yum. The ionized gas traced by Bry distributes 
along ~7'.'5 of slit #1 and ~4" of slit #2. Similarly to flie opti- 
cal observations, the two-dimensional NIR spectrum along slit 
#1 shows two peaks in the surface brightness distribution of 
the emission lines that correspond to the NE blob and the SW 
arc. Also within these peaks, we detect a featureless and weak 
continuum, specifically at the (0",0") and (0",-5") locations. 
Spectra extracted from 0'.'75x0'.'75 regions around these con- 
tinua are shown in Fig. [3] 

We investigated the nature of the continua detected in 
the /T-band. NIR continuum can be produced not only by a 
stellar photosphere but also by: (1) the free-free and free- 
bound emission from ionized gas, (2) hot dust and (3) scat- 
tered light, i.e. radiation escaping through optically thin paths 
of the nebula and reflected into the line of sight of the ob- 
server The contribution by free-free emission to the observed 
infrared flux density can be est imated using equation 7 of 
iRom an-Lopes & Abraham ( 20041 ) and the free-free Gaunt fac- 
tors of .Hummer (,1988i) . For a flux density of 884 mJy at 5 GHz 







Line fluxes (10 ' 


** W m-2) 


A \J1IU) 


Identincation 


out ff I 




2.033 


H2 1-0 S(2) 


28.3 + 3.2 


3.2 ± 1.0 


2.042 


He 6^P-43S 


7.7 + 1.5 


2.3 ± 0.5 


2.058 


He2'P-2'S 


1406 ± 39 


341 ± 10 


2.073 


H2 2-1 S(3) 


15.0 ± 2.4 


< 3 


2.112 


He 43 'S-3' 'P 


84.7 ± 4.0 


20.9 ± 1.2 


2.121 


H2 1-0 S(l) 


147.1 ± 3.5 


20.5 ± 0.9 


2.145 


[Fein] ^G3-^H4 


15.2 + 1.4 


3.5 ± 0.6 


2.154 


H2 2-1 S(2) 


<7 


< 2 


2.161 


He73'F-43'D 


61.1 ± 3.3 


16.0 ± 1.1 


2.165 


Bry 


1994 ± 106 


495 ± 27 


2.198 


UIR 1 


13.1 ±2.5 


2.6 ± 0.5 


2.218 


[Fern] 'Gs-'He 


63.4 + 1.7 


15.7 ± 0.8 


2.223 


H2 1-0 S(0) 


45.6 ± 2.2 


7.1 ± 1.0 


2.242 


[Fern] ^G4-'H4 


25.5 ± 1.3 


7.6 ± 1.3 


2.247 


H2 2-1 S(l) 


14.3 ± 1.1 


< 5 


2.286 


H2 3-2 S(2) 


15.8 ± 1.8 


4.0 ± 1.0 


2.340 


Hi 31-5 


11.4 ±3.4 


3.9 ± 1.2 


2.344 


Hi 30-5 


16.0 ± 4.8 


< 8 


2.349 


Hi 29-5 


60.7 ± 18.2 


12.5 ± 3.7 


2.354 


Hi 28-5 


16.4 ± 4.9 


3.1 ±0.9 


2.360 


Hi 27-5 


20.1 ± 6.0 


3.8 ± 1.1 


2.367 


Hi 26-5 


28.0 ± 8.4 


5.1 ± 1.5 


2.374 


Hi 25-5 


19.7 ±5.9 


5.0 ± 1.5 


2.383 


Hi 24-5 


51.4 ± 15.4 


10.4 ± 3.1 



Errors are at the Icr level and only include the uncertainties coming 
from the fitting procedure. In the case of the H i Pfund series, we give 
a larger error of the order of 30% because of the uncertainty fitting 
the continuum. The H? 3-2 S(3) transition at 2.201 ^m is observed 
in regions 1 and 4 of slit #1 (see Table O. This line is probably 
blended with a transition of [Fe m]. This line is likely contaminated 
by atmospheric lines. 



(IZiilstra et al. . I1990I) and an electron temperature of 10000 K, 
the expected flux density in the /T-band is about 130 Jy, which 
translates into ~ 8 x 10"'"* Wm"^/im"' at 2.2 fim. This is com- 
parable to the continuum flux level detected across slit#l (cf. 
Fig.|3]l. It seems, therefore, that free-fee emission does indeed 
dominate the continuum observed in the NIR and hides the stel- 
lar photospheres detected in the optical within the NE blob and 
the SW arc. 

Moving along slit #1 and from NE to SW, we only start 
detecting emission from Bry at the (0",0") position, i.e. the lo- 
cation of the continuum within the NE blob. After this position, 
the Bry line flux slowly increases, reaching a local maximum 
r.'25 ahead of the NE blob and coincident with the weaker blob 
detected in the acquisition image. It then reaches a minimum 
within the nebula and increases rapidly afterwards. This bright 
portion corresponds to the SW arc structure seen in the image 
and reaches an absolute maximum 5'.'5 ahead of the NE blob. 

Moving along slit #2, the Bry emission resembles that of 
a shell structure, with two maxima at the edge of the ionized 
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Fig. 3. (Left) Spectra of Ml-78 extracted across (top) slit #1 between -7'.'5 and 0" and (bottom) slit #2 between -2" and 2"(see 
Fig.[T]). The observed lines, their identification and their line fluxes are listed in Table[3] (Right) Spectra extracted across slit #1 
from (y.'75x(y.'75 regions around the continua detected within the NE blob (top) and the SW arc (bottom). Note the logarithmic 
scale of the flux axis. 



nebula (at about - 
the nebula. 



1" and 2'.'5) and a decreased emission within 



1.0 



The following panels compare the He i distributions of the 
2.058 and 2.112 jim lines to that of the Bry. It is noticeable that 
while the 2.112 fj.m Hei line spatial variation traces the Bry 
emission well, this is not the case of the 2.058 fim Hei line. 
This difference will be discussed in detail in Sect. 14.41 



The spatial variation in the H2 1-0 S(l) transition at 
2.121 fim is also plotted in Fig. |4] The H2 emission extends 
over ~30" across slit #1 and ~20"across slit #2. We distin- 
guised four different regions in slit #1 and three in slit #2, which 
will be analysed separately in Sect. 14.71 

Ml-78 has also been observed using the mid-infrared long- 
slit spectrograph SpectroCam-10 mou nted on the 5 m Hale 
Telescope at Palomar Observatory (cf. iPeeters , I2OO2'). These 
observations were taken using a slit of l"xl6'.'4 slit positioned 
along the same axis of our slit #1 and centered on the brightest 
spot. Figure|5]compares the spatial variations of the [Neii] line 
at 12.8 fim and the flux at 12 fim with the Bry profile. Both sil- 
icate in emission and dust continuum contribute to the flux at 
12 fim. The variation of the [Ne 11] line traces that of Bry well, 
with only a small variation close to the NE blob. The "dust" 
profile shows a cavity within the nebula and a clear shell-like 
structure, suggesting that the dust have been removed from 
within the ionized gas volume either through radiation pressure 
from the star within the NE blob or destruction. In contrast, in 
the arc, the "dust" emission concides with that of the ionized 
gas. 
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Fig. 5. Spatial variation of the [Nen] 12.8 fim fine-structure 
line (solid line) and the 12 fim dust continuum (dashed line) 
obtained with the mid-infrared spectrograph SpectroCam-10 
mounted on the 5 m Hale Telescope at Palomar Observatory. 
These variations correspond to a l"xl6'.'4 slit positioned along 
the same axis of our slit #1. We show the Bry variation (dotted 
line) for comparison. 



4.3. Unidentified features 

Several unidentified features are also present. Two relatively 
common unidentified features at 2.1981 and 2.2860 ^m are 
present in the K-han d spectra of e.g. O r ion ( cf. Luhman et al. , 
[l998i) and PNe (cf. iLuhman & RiekS. [l99fe ,Hora & Latter . 
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Fig. 4. (Left) Spatial variation of Bry and the two strongest He: lines at 2.058 and 2.112 /im across slit #1 (left) and slit #2 
(right). The normalized He i intensity distributions are compared to the Bry profile (dashed line). (Right) Spatial variation of 
the H2 1-0 S(l) transition at 2.121 /im across slit #1 (top) and slit #2 (bottom) compared to the Bry profile (dashed line). The 
different regions analysed in the text are identified (dotted lines). The origin corresponds to the (0",0") position of Fig. [T] 



19961: iHora et al.1 Il999t iLumsden et aU bOOl b"). quite close 
to the 3-2 S(3) and 3-2 S(2) transitions of H2 at 2.2008 and 
2.2864 fim, respectively, wher e these are the air wa velengths. 
This is discussed at length by Geba lle et al.l (Il99lb . who ex- 



amined lists of permitted and forbidden atomic lines for wave- 
length coincidences with these features and only found near 
coincidences with weak lines. Geballe et al. claimed that the 
parent ion responsible for these features corresponds to an ex- 
citation level between 30 and 60 eV, but the higher end of 
this range is inconsistent with the observed presence of these 
lines in objects with e ffective temperatu res < 40000 K (cf. 
Lumsden et al.l 200 1 b ) . [Pinersteiri (12001 ) claimed an identifi- 
cation for these features as [Krm] and [Seiv] for the 2.198 and 
2.286 fim lines, respectively. However, the abundance of both 
krypton and selenium in the interstellar medium is sufficiently 
low that this remains a mystery. 

In the integrated spectra shown in Fig. [3] we detect two 
Unes at 2.1981 + 0.0001 yum and 2.2864 + 0.0001 yum. While 
the former line is coincident with the unidentified feature dis- 
cussed by Geballe et al. (and designated UIR 1), the other line 
can be unquestionably identified as the 3-2 S(2) transition of 
H2. We however note that the 3-2 S(3) transition of H2 (with 
an observed central wavelength at 2.2008 + 0.0001 yum) is de- 



tected in regions 1 and 4 of slit #1, well outside the ionized 
gas region (cf. Fig. HJ. In these two regions, the UIR 1 feature 
is not detected. The UIR 1 feature is not detected either in re- 
gions 1 and 3 of slit #2. This behavior suggests that the carrier 
of this unidentified UIR 1 line is probably an ion created by UV 
photons of at least 13.6 eV. 



4.4. Analysis of the He 1 recombination lines 

Calculations of the He i recombination ca scade spectrum have 
beer i published by differen t autho rs, e.g. Smits ( 1991 . 19961) 
and Beniamin et al. (19991 l2002h . Each study improved the 



prior treatment of physical processes, mainly as the result of 
improved theoretical calculations of various rates. The He i re- 
combination cascade spectrum has been recently revisited by 
Bauman et al. (2005) and Porter et al.. (2005, 2007) using im- 
proved radiative and colUsional data. They use highly accurate 
calculations of the 7-resolved transition probabilities that take 
singlet-triplet mixing into account. Calculations of collision- 
les s (low-density l imit) C ase B He: emissivities are published 



by iBauman et al. I tP05). The source of the 7-resolved code 



used in these calculations is available online (see iPorter..20()7l) . 
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Theoretical predictions of several He i line r a tios ar e shown 
in Fig. |6] We used the calculations of iPorteii (120071) . which 
are only avai lable for the low densit y limit and for T3889 - 0, 
and those of iBeniamin et al. whose predictions con- 



sider different values of Hg and T^g. Here, Tjggg is the line cen- 
ter optical depth of the 3^P-2^S 3890 A line. Since the ratios 
are rather insensitive t o tig, we only represented the curves of 
Benjamin et al. ( 1999h for ng — 10'^ cm"^ but varying from 
5000 to 12500 K and 73889 from to 100. Fr om the figure 
it can be seen that the theoretical predictions of lPorteii (120071) 
with Je between 8000 K an d 10000 K reproduce the data well. 



whereas the calculations of Benjamin et al. d 19991) give temper 



atures somewhat lower than 8000 K for T3889 = 10. These val- 
ues of Tf. - corresponding to that characteristic of the He^ zone 
of the nebula - are lower than those obtained from the [N 11] 
line ratios. This is a rather common behavior in Galactic HII 



regio ns, as can be found in the Orion Nebula dEsteban et al 



20041) . NGC 3576 dOarcfa-Rojas et aU 120041) and others. The 
difi^erence between Te obtained from coUisionally excited and 
recombination lines could be related to th e presence of tem per- 
ature fluctuations in ionized nebulae (e.g. lPeimbert . 1967 ). 

Figure |4] shows that while the He: 2.112/Br7 ratio is 
fairly constant across the slits, there nevertheless exists a large 
spatial variation in the Hei 2.058//m line with respect to 
Bry and the other He i recombination lines. For instance, the 
He 1 2.058/2. 112 ratio varies between 10 and 25, reaching its 
lowest values at positions < -5"for slit #1 and > O"for sUt 
#2. The origin of this spatial variation is not clear. To first or- 
der, temperature variations should have similar effects on all 
the lines, so this is probably not the cause. Density variations, 
on the other hand, might be able to produce the observed effect 
but would in principle not be measureable. The most probably 
cause of this effect seems to be optical depth variations across 
the nebula. 



The theoretical calculations of iBauman et al.l (l2005h fail 
to reproduce the He: 2.058//m line. For instance, they pre- 
dict a He 1 2.058/2. 1 12 ratio around 0.3-0.4. Evidently, the line 
flux of th e He I 2.058 Aim line is largely underpredicted. The 
models of iBenjamin et al. (Il999h . on the other hand, predict a 
He 1 2.058/2. 1 12 ratio between 20 and 50 for large variations of 
T3889, density and temperature. Hence, models seem incapable 
of fitting the He 1 2.058 /urn line. This is because this line is not 
well defined in the Case B approximation: e.g. the Bauman et 
al. code assumes that the 2'P to I'S photons are absorbed by 
some other atom or ion. The He: 2.058 jum transition is de- 
termined by the population in the 2'P state. The Case B ap- 
proximation assumes that 2'P to 1 'S photons do not escape the 
cloud, but does that mean that they are reabsorbed and com- 
pletely converted to 2'P to 2'S photons, or that they are de- 
stroyed or absorbed by some other atom or ion? This effect, 
negligible for every other line, has an enormous effect on the 
2.058 jum line. 

4.5. The H 1 Pfund lines 

The forest of lines redward of about 2.3 ;um corresponds to 
the higher members of the atomic Hi Pfund series. We de- 



tected H I pfund series lines from Pf24 to Pf3 1 which allowed 
us to infer the extinction across the A"-window when compared 
with the theoretical pr edictions of iHummer & Storey! ( 119871) 
and I S tore V & Hummer! (11995). These lines are consistent with 
an extinction A/f < 1 in agreement with the optical estimate. 

4.6. Forbidden iron lines 

We report the detection of [Fe iii] emission. The ^T-band hosts 
[Fe III] lines from the multiplet connecting the and ^H terms 
(e.g. figure 6 of Bautista & Pradhan, 1998). Of the four transi- 
tions in this multiplet with relatively large (magnetic dipole) 
transition probabilities, three are in this wavelength region: 
2.145, 2.218 and 2.242 fim. The fluxes measured for these lines 
are fisted in Table|3] Since these lines originate from highly ex- 
cited levels, their absolute emissivities are quite sensitive to the 
temperature. However, the fine ratios are primarily sensitive to 
the density. 

We compared our observed [Feiii] line rati os with theo 



retical predictions. We used the calculations of iKeenan et al 



( 2001!) and the mo re complete ones provided by M. Rodriguez 
(Rod riguezl!2002l) . These last calculations use a 33-level model 
atom where all collisio n strengths are those calculated by 
Zhang & Pradhan j 19971). the transition probabilities those rec- 
ommended bylfioese Fischer & Rubliil ( 1998 ) (and those from 
Garstangl 19581 for the transitions not considered by Froese 
Fischer & Rubin), and the level energies have been taken from 
the NIST database. We considered densities between 10^ and 
10^ cm"^. In Fig. |2l we can see that our data points are con- 
sistent with the predictions of Rodriguez, giving densities on 
the order of those determined from the optical forbidden lines. 
However, the predictions of Keenan et al. fail to fit the observed 
2.242/2.218 ti m ratios; this is also the ca se for the theoretical 
predictions of Bautista & Pradhan (!l998'). who overpredict the 
observed 2.242/2.218 jum ratios by factors between ~1.5 and 
2.5. 



4.7. H2 excitation 

The NIR emission from H2 can be produced by either ther- 
mal emission in shock fronts or fluorescence excitation by 
non-ionizing UV photons in the Lyman- Werner band (912- 
1108 A). These mechanisms can be distinguished since they 
preferentially populate different levels producing different 
spectra. For shock-excited H2 the lower energy levels are typi- 
cally populated as for a gas in local thermal equilibrium (LTE) 
with a characterisitic excitation temperature around lOOOK. 
In contrast, for radiatively excited gas the population follow a 
non-LTE distribution characterized by high excitation temper- 
atures around lOOOOK and lines from high vibrational levels 
(v » 1) may be detected. 

Non-thermal excitation mechanisms readily excite the v - 
2 and higher vibrational states, whereby collisional transi- 
tions preferentially de-excite the v - 2 level in favor of 
the V = 1 level. Hence, the H2 1-0 S(l)/2-l S(l) ratio 
has on occasion been used to distinguish between these two 
processes. However, in dense clouds, FUV irradiation can 
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Table 4. Air wavelengths, upper energy levels, Eistein A coefficients, degeneracies and line fluxes of the H2 lines detected in 
Ml-78. We also give significant line flux ratios and results such as the ortho-to-para ratio (0) and the excitation temperature. 



Transition 








81 






Slit 1 






Slit 2 






(jim) 


(K) (10-^ s) 




Region 1 


Rej 




T? cirri Oil ^ 


Region 4 


Region 1 






Region 3 


1-0 S(2) 


2.0332 


7584 


3.98 


9 


35.3 ± 2.6 


18.5 


± 1.9 


10.4 ± 1.1 


26.3 ± 1.1 


8.0 ± 1.4 


2.9 


±0.7 


5.3 ± 1.0 


2-1 S(3) 


2.0729 


13890 


5.77 


33 


16.7 ± 1.9 


9.9 


± 1.1 


4.6 ± 0.6 


22.3 ± 1.1 


5.9 ± 1.1 




< 3 


3.1 ± 0.7 


1-0 S(l) 


2.1212 




^ 47 


21 


148.6 + 10.4 


92.3 


± 1.4 


48.5 ± 1.2 


150.8 ± 3.6 


44.5 ± 2.3 


19.1 


± 0.8 


30.1 ± 1.9 


2-1 S(2) 


2.1536 


13152 


5.60 


9 


< 3.5 




<4.2 


< 2.5 


10.7 ± 1.1 


3.6 ± 0.9 




< 2 


< 3 


3-2 S(3) 


2.2008 


19086 


5.63 


33 


5.7 ± 1.4 




< 4.2 


< 2.2 


8.5 ± 2.3 


< 3 




< 3 


< 2 


1-0 S(0) 


2.2229 


6471 


2.53 


5 


33.8 ±3.8 


28.5 


±2.2 


16.6 ±0.8 


38.3 ± 1.7 


11.6 ± 1.5 


6.9 


± 1.0 


< 5 


2-1 S(l) 


2.2471 


12550 


4.98 


21 


15.1 ± 2.2 


9.3 


± 1.1 


4.2 ± 0.7 


16.3 ± 1.6 


5.0 ± 1.1 




< 4 


4.0 ± 1.1 


3-2 S(2) 


2.2864 


18386 


5.63 


9 


< 3.3 


11.6 


± 1.3 


4.0 ± 0.6 


<4.5 


< 3 


4.4 


± 1.0 


2.3 ± 1.0 


1-0 S(l)/2-l S(l) 








9.8 ± 1.6 


9.9 


± 1.2 


11.5 ± 1.9 


9.2 ± 0.9 


8.6 ± 1.9 




> 4 


7.5 ± 2.1 


1-0 S(l)/3-2S(3) 








26.1 ± 6.6 




> 22 


> 22 


17.6 ±4.8 


> 17 




> 7 


> 15 


01 










3.5 ± 0.3 


3.1 


± 0.2 


2.9 ± 0.2 


3.7 ±0.1 


3.6 ±0.4 


3.3 


± 0.4 


(b) 


02 










A 




A 


A 


2.0 ± 0.2 


1.7 ±0.4 




(f) 


(<-■) 


(K) 










980 ± 165 


720 


±55 


690 ± 20 


720 ± 150 


750 ± 175 


570 


±55 


590 ± 120 


Tex (K) 










2255 ± 230 


2100 


t260 


1910 ±330 


2465 ± 390 


2380 ± 330 




id) 


2265 ± 420 



Line fluxes are given in units of 10""* Wm"^. The 1-0 S(0) line was not been measured. Only one or two v = 2 lines were detected. 
There are only upper limits for the v = 2 transitions. Excitation temperature determined from the v = 1 states only. Excitation temperature 
determined from the v = 1 and v = 2 states. 
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= 100 




- 5000 K A ' 
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Fig.6. Plot of the Hei 2.1 1243-'S-3^'ip/2.161 T^'F-^^-'D 
ratio against the He: 2.0426^P-4^S/2.11243'iS-3^''P ratio. 
Ratios measured for slit #1 are plotted by a black square, 
and by a black triangle for slit #2. The open triangles repre- 
sent the theoretical p redictions using the /-resolved code of 
iBauman et al.l (|2005|) for = 5000 K, 8000 K, 10000 K and 
12 500 K. Similarly, the o pen circles represent the predictions 
of iBeniamin et al. (1 19991) for the same electron temperatures 
and three values of T3889: 0, 10 and 100. 



heat the gas to temperatures ~1000K leading to collision- 
ally excited and thermalized 1-0 line emission, while direct 
FUV pumping maintains non-thermal fluorescent populations 
in higher-lying vibrational lev els (e.g. Sternberg & Dalgarnol 
1989l:lDraine & BertoldiL ll996V Hence, for low densities (hh ~ 




0.2 0.4 

[Fe in] 2.145/2.218 



0.6 



Fig. 7. .g-b and [Feiiil lin e ratios . We show the theo retical pre- 
diction s of Keenan et al. (2001) (dotted line) and Rodrfgue^ 
(I2OO2I) (dashed line) for densities between logWe - 2 and 
log^e - 5. The black square and triangle represent the line 
ratios measured for slit #1 and slit #2, respectively. 



10^ cm ^), the 1-0 S(l)/2-l S(l) ratio retains a pure fluores- 



cence value around 2. But for high density clouds (nu > 10^ - 
10"^ cm"^) and an intense incident FUV radiation field with 
Go ^ 10"* (where Go is the incident FUV photon flux between 
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6 and 13.6 eV measured in units of the local interstellar radi- 
ation field), the 1-0 S(l)/2-l S(l ) ratio will rapidily approac h 



_ap 

a thermal value of about 10 (e.g. iHoUenbach & Nattalll995l) . 



which has often been interpreted to indicate shock emission. 
Conversely, the 1-0 S(l)/3-2 S(3) ratio retains a value of ap- 
proximately 8 in a dense photodissociation region (PDR) when 
Go ~ 10^ and only for very high densities and a very intense 
FUV radiation field («h > 10^ cm^^ and Go > 10"^) will this ra- 
tio a pproach a thermal value of 10-100 ( Hollenbach & NattaL 
1 19951) . Therefore, the combination of the 1-0 S(l)/2-l S(l) and 
1-0 S(l)/3-2 S(3) ratios should in principle serve as an indica- 
tor of the excitation mechanism of the H2. 

We distinguished four different regions along slit #1 and 
three along slit #2 with sufficient signal-to-noise (cf. Fig. HJi. 
The line fluxes measured for all these regions are listed in 
Tableg] This table also lists the values for the 1-0 S(l)/2-l S(l) 
and 1-0 S(l)/3-2 S(3) ratios. Across slit #1, the 1-0 S(l)/2-l 
S(l) remains practically constant within the errors with an av- 
erage value around 9.7 + 0.6. The 1-0 S(l)/3-2 S(3) ratio is in 
general higher than 20 for all these regions. A similar scenario 
occurs for the three regions along slit #2. Both ratios suggest 
thermal emission. However, in the case of a dense PDR with 
«H ^ 10^ cm"^ and a very intense FUV radiation field with 
Go > lO"*, these ratios can be explained by a non-thermal exci- 
tation mechanism. We can now estimate if the PDR in Ml-78 
can effectively be characterized by such high values of Go and 

"H- 

Go can be derived by assuming that all the UV pho- 
tons are absorbed in a spherical shell of the size of 
the PDR and re-emitted in the infrared continuum (cf. 
Bernard-Salas & TielensL 2005 ). Normalizing Go to the aver- 
age interstellar field (1.6 x 10"^ W m"^; .Habing. .1968.) . the 
expression is: 



4F, 



IR 



1.6 X 10-6 X 2.35 xlO-i'6l2 



(1) 



where is the diameter of the PDR in arcseconds and Fir is the 
observed infrared flux in W m"^. We use a FIR flux of ~ 7 x 
10 W m"^ determined by fitting a blackbody t hroug h the 



IRAS fluxes at 25, 60 and 100 fim (iPeeters et al 
obtain a Go around (1 - 2) x 10^ considering a size of 20 



20021). We 

30" 



(see Fig.m. We note, however, that Go might be substantially 
larger if the filling factor of the PDR gas is small. 

The total density can be determined, for instance, 
using fine-structure lines emitted in the PDR. The Oi 
63//m/[Cii]158yum ratio is a good tracer for densities in be- 
tween lO"' and 10* cm"^ because these two lines have different 
critical densities. These two lines have been mea sured by ISO 
within an aperture of ~ 80" /Peeters et al.','2002'). With a ratio 
of 7.5 + 0.6, the models by Kaufman et al^ (1999.) indicate a 
density of 10"* - 10^ cm"^ when Go > lOl 

Therefore, it seems reasonable to assume the case of a very 
dense PDR and a strong radiation field. This analysis, however, 
is not conclusive and we will use two additional diagnostics to 
clarify the origin of the H2 emission: the ortho-to-para ratio and 
the excitation diagrams. 



4.7.1 . The ortho-to-para ratio 

An additional diagnostic is provided by the ortho-to-para ra- 
tio. At the moment of formation, molecular hydrogen can join 
one of two denominations: ortho (aligned nuclear spins) or para 
(opposed nuclear spins). The ortho-to-para ratio (p, based on 
the statistical weights of the nuclear spins, is 3:1. In outflow 
regions where shock excitation is the primary emission mecha- 
nism, observations of vibrational emission lines typically re- 
veal ortho-to-para ratios f or vibrationally ex cited states that 



are comparable to 3 (e.g. ISmith et al. . Il997l) . as is expected 



if the gas behind such shocks is in LTE. Ho wever, as dis- 
cussed at length by Sternberg & Neufeld ( 19991) . the ortho-to- 
para ratio can take values other than the intrinsic 3 if the exci- 
tation is caused by UV fluorescence, since the transitions are 
pumped by optically thick UV transitions in which the optical 
depth depends on whether the molecule is in an ortho or para 
state. Indeed, ratios in the range 1 .5-2.2 have been measured 
in PDRs (e.g. [Ra msay et al., 1993: IChrvsostomou et al. . 1993 
Shupe et al.Ul998tlLumsden et alll2001ah . 



From the 1-0 S(0), 1-0 S(l) and 1-0 S(2) line fluxes, we 
measured the ortho-t o-para ratio for the v = 1 states (0i) fol- 
lowing equation 9 of ISmith et alJ (119971) . In Table|4] we list the 
V = 1 ortho-to-para ratio for each region. These are consistent 
with a ratio of 3 suggesting a thermalization of the v = 1 states. 
Of course, they do not tell the whole story. In a few cases, a 
value of the ortho-to-para ratio for the v - 2 states (^2) was de- 
rived from the 2-1 S(l), 2-1 S(2) and 2-1 S(3) transitions using 
equation 1 of Davis et al.l (l2003h . We obtained ratios around 2, 
typical of PDRs and in agreement with the ratios derived by 
Lumsden etall (l2001bl) . 



4.7.2. H2 excitation diagrams 

The simplest way of characterizing the molecular hydrogen 
emission, however, is to plot the observed column density 
against the energy of the upper level. The measured intensity, 
/, of a given H2 line can be used to calculate the column den- 
sity of the upper excitation level of the transition, which, for 
optically thin emission, is given by: 



Ni 



AnA£ 
Ajhc 



(2) 



where Aj is the rest wavelength and A, is the Einstein A- 
coefficient taken from Turner etal. (1977). If collisional de- 
excitation is assumed to dominate, the H2 will be in LTE and 
the energy levels will be populated in a Boltzmann distribution. 
The relative column densities of any two excitation levels can 
thus be expressed in terms of an excitation temperature Tex: 



— - — exp 



Sj 



-(El - Ej) 



(3) 



where gj is the degeneracy, Ej is the energy of the upper level 
taken from iDabrowskil (119841) and k is the Boltzmann's con- 



stant. The values of Aj, Ej, Aj and gj for the lines detected in the 
spectra are shown in Table ID Clearly, plotting the logarithm of 
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Fig. 8. H2 excitation diagrams for each region. The lines represents the best fitting purely thermal single-temperature models. 
The solid line is the fit to the v = 1 level only and the dotted line is the fit to all the data except the v = 3 level. 



the ratio of the column densities divided by the appropiate sta- 
tistical weights gives the excitation temperature from the slope 
of the data. 

In Fig. [8] we show plots of logNj/gj versus the energy 
for each region. We normalized the population distributions 
relative to that inferred from the 1-0 S(l) line. For ther- 
malized populations at a fixed gas temperature, the logN/g 
points should lie on a straight line in these excitation dia- 
grams. Cleary, such single-component models do not fit the 
data. While the v = 1 levels appear to be thermalized at 
r ~ 600 - lOOOK (solid fines), the inclusion of the v = 2 
levels raises t he best-fitting te r nperatu re toT ~ 2000 - 2500 K 
(dotted lines). iLumsden et al. I (l2001bh find excitation tempera- 
tures of 2400 K at the position of the nebula and slightly higher 
(around 3500 K) to the north and south. 

However, even at these high temperatures, the v = 3 levels 
remain underpopulated and the excitation temperature inferred 
from the relative populations in the v = 2 and v = 3 levels ex- 
ceeds 4000 K. At such gas temperatures, the molecules would 
be rapidly dissociated, suggesting that a non-thermal excitation 
mechanims is responsible for the excitation of the v = 2 and 
V = 3 levels. 

The results of these diagrams indicate again that the ob- 
served H2 excitation is produced in a dense PDR rather than in 
shocks. In fact, a high vibrational excitation temperature, such 
as the 4000 K inferred from the relative v - 2 and v = 3 pop- 
ulations, is a basic characteristic of FUV-pumped fluorescent 
emissi on (e.g.lBurton et al. [ Il990l) . Furthermore, and as pointec 
out by Davies et alJ ( 2003 ). an additional signature of fluores- 
cent excitation appears to be present in the data. For those re- 



gions in which the 2-1 S(l), S(2) and S(3) lines are detected 
(region 4 in slit #1 and region 1 in slit #2), the resulting N/g val- 
ues for the j-3, 4 and 5 levels do not appear to be thermalized. 
Instead, the y = 4 para-level is shifted above the adjacent j-3 
and 7 = 5 ortho-levels. The ortho-to-para ratio for the v = 2 
transitions appears to be out of equilibrium and suppressed be- 
low the value of 3 obtained in LTE. This is precisely what is 
expected in the case of FUV pumping, where the populations 
of the vibrationally p umped ortho-H^ is suppress ed relative to 
those of the para-H2 ( [Sternberg & Neufeld , Indeed, we 

find values as low as 2 for the v = 2 ortho-to-para ratio. 



5. Discussion: tlie nature of l\/l1-78 



Ziilsfi-a et all (Il990l) and IPerinottol (Il99lb reported that the 



chemical composition of Ml-78 is enriched in nitrogen with 
respect to oxygen when compared with the typical abundance 
pattern of Hii regions. Our optical observations, much deeper 
than previous ones, confirm this enrichment but also indicate 
that it seems to be stronger in the NE blob. The blob coin- 
cides with a very faint stellar emission that can be seen in our 
red optical (as discussed in Sect. 14.21 the NIR continuum is 
on the other hand dominated by free-free nebular emission). 
Therefore, it seems that a process of localized chemical en- 
richment is taking place to the NE of Ml-78, and that the star 
we see at the blob is its most probable source. At the blob, 
the nitrogen enrichment is correlated with a defficiency in the 
O abundance and a - dubiou - He enrichment. This abundance 
pattern is typical of nebular ejecta around evolved massive stars 
such as WR (e.g. .Esteban et"all 1992) and LBV stars (e.g. 
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Lamers et al.Ll200ll) . The chemical composition of these nebu- 



lae reflects the product of the CNO nucleosynthesis in the H- 
burning zones of the parent post-main-sequence massive star. 

In order to explore the likeliness of stellar ejecta in Ml- 
78, we derived its ioni zed mass from the ra dio flux of 884 mJy 
at 5 GHZ obtained bv lZiilstra et al. Considering a dis- 

tance of 8.9 kpc and a diameter of 8'.'0, we obtain ne(rms) ~ 
6200 cm"^ and M(Hii) ~ 3 Mq. The fact that the ne(rms) coin- 
cides with the «e obtained from the optical forbidden lines indi- 
cates that the filling factor of the gas is very close to 1 . The ion- 
ized mass derived for Ml-78 is on the order of those found for 
ejecta nebulae arou nd WR and LBV star s , typically about a few 



solar masses (e.g. lEsteban et all Il992t iLamers et all 1200 lb 



However, our nominal value of 3 Mg should be considered as 
an upper limit of the mass of the ejecta. 

The spatial variations of the physical conditions and chemi- 
cal abundances and the presence of more than one possible ion- 
izing source indicate that Ml-78 is better described as a combi- 
nation of a compact H ii region + ejecta. In this se nse, this com- 
bined nature can explain the negative result of iHutsemekers 
(Il997h . who doubts the classification of Ml-78 as an ejecta 
nebula due to its very large dust mass obtained from the far- 
infrared emission measured by the I RAS satellite. We ha ve to 
consider that the ejecta nebulae that HutsemekersI ( 1997 ) uses 
for his comparisons are isolated and are not related to a nearby 
star-forming region as Ml-78 presumably is. Further indirect 
evidence of an ejecta nebula in Ml-78 comes from the afore- 
mentio ned detectio n of strong [Ni ii] emission in its spectrum. 
In fact. lLucvl(ll995h indicates that circumstellar nebulae around 
LBV stars are anomalous [ Nin] emitters, as has b een found in 



the nebula around P Cvgni ([Johnson et al.L 119921) and two ob 



jects in flie LMC JStahl & Wold.ll986h . However, die high ex- 
citation and of Ml-78 do not support its identification as a 
LBV n ebula. Typically, thes e 



K (e.g. iLamers et al.L 12001 



or even unobservable [O iii] (iNota et al 



objects sho w lower than 10000 
1995), and very faint 
19951) . Obviously, this 



Nota et al. 



difficulty can be obviated if there is another hotter ionizing 
source in the nebula apart of the hypothetical LBV star itself. 
In contrast, ring nebulae around WR stars show a wide range 
of nebular excitations and Tg, basically depending on the spec- 
tral classification of the ionizing WR star. In fact, between 
12000 and 14000 K have been found in nitrogen-rich nebulae 
ionized by early WN and WO stars as NGC 2359, S 308 or 
G 2.4-1-1.4 dEstebanet al.lll990lll992h . 

Another remarkable feature of Ml-78 is the expansion ve- 
locity of the ionized g as of about 25 km s"' measured by 
Gussie & Tavloii (119891) . Although this fact was considered by 
Gussiel(ll995 ir as an indication of the PN nature of the object. 



such supersonic e xpansion velocities are also common in WR 



ring nebulae (e .g. Rosadol 19861) as well as in LBV nebulae 



fe.g. lNota et al 



1995). 



A close analogue of 



(IHeydari-Malayeri & Testoil Il982h . 



Ml-78 
a compact 



IS 

Hii 



N82 
region 



in the Large Magellanic Cloud. This object shows nitrogen 
enrichment and oxygen defficiency. The ejecta seem to be 
produced by a WC 9 star and ther e is one additional ionizing O 
star in the nebula /Moffat' 1991, however, argued against this 
classification and suggested that the presence of a WNL or Of 



star in N82 is more likely). The linear size, ionized mass and 
He of N82 (0.4 pc, 7 Mq and 5500 cm"-*, respectively) are very 
similar to those of Ml-78 but the higher excitation and of 
this Galactic object would indicate that the progenitor should 
be a WC or WN star of an earlier - hotter - subclass. 

We investigate now the possibility that a second star, specif- 
ically an O star, might be present in Ml-78, e.g. at the position 
of the SW arc where a second stellar continuum is detected 
in the optical. An estimate of the spectral type of the ionizing 
star in H ii regions based on the He i 2. 1 12/Bry ratio has been 
proposed by Hanson et al. (2002). This method relies on a mea- 
sure of the ionization fraction of helium in the H n region. For 
very hot stars (Teff > 40 000 K), helium is ionized throughout 
the entire nebula producing a constant relative strength of the 
recombination lines of helium and hydrogen. For a Teff below 
this value, the helium emission compared to hydrogen quickly 
decreases as the H ii region region be comes predominantly neu- 
tral helium. Using the predictions of iBauman et alJ (l2005l) . we 
determined the ratio of He i 2. 1 12 fim to Bry when the helium 
is fully ionized in an Hii region. Assuming «(He)/n(H) - 0.10 
and an electron temperature of lOOOOK, these models predict 
a ratio of 0.040, in perfect agreement with the observed ratio in 
both slits (0.042 + 0.002, see Table©. Assuming for now that 
a single O star is responsible for the ionization of the nebula, 
this would imply that the ionizing star must be an early O star 
with Teff > 40 000 K. This would correspond to a spectral type 
earlier than 05.5 V, as inferred from the new calibration of O 



star parameters published by' Martins et al.l (l2005h and based on 
state-of-the-art stellar models. 

We do indeed find a He 1 2. 1 12/Bry ratio of 0.043 + 0.003 
at the position of the SW arc. However, this ratio is lower 
at the position of the NE blob; specifically, Hei2.112/Br7 = 
0.029 + 0.005 in the 0'.'75x0'.'75 region around the steUar con- 
tinuum within the NE blob. This implies either a low electron 
temperature (< 5000 K) in the NE blob (which is not the case, 
see Section[32]i or the presence of a colder star that is not able 
to fully ionize He. It therefore appears that the WR star that we 
suspect to be hiding within the NE blob is not solely responsi- 
ble for the ionization of the nebula. 

It seems reasonable, then, to assume that a early O star is 
present in the SW arc of Ml-78. The NIR spectrum of such 
a massive star will be characterized by emission of [N iii] at 
2.115 jum but the detection of this line at intermediate spec- 
tral resolutions is difficult since it blends with the He i line at 
2.112 fim. Moreover, it is not surprising that we have not de- 
tected direct signs of its presence, since the NIR emission we 
detect is merely nebular Detection of stellar features in the op- 
tical spectrum is also difficult because of the high visual extinc- 
tion. 

The spectral type of the ionizing stars in Hii regions 
can be also determi ned from radio continuum measurements. 
Zijlstra et al. d 19901) measured 884 mJy at 5 GHz. Assuming 
that the nebula is dust-free, optically thin, homogeneous and 
ionization-bounded, this density flux gives a minimum value 
for logQo of 48.80 dex, where Qo is the rate of Lyman contin- 
uum photons required to maintain the ionization of the nebula. 
This value of go would correspond to a star with a spectral type 
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earlier than 06.5 V star, in perfect agreement with the previous 
estimate. 



6. Conclusions 

There is considerable controversy surrounding the nature of 
Ml-78, a compact object located beyond both the Local arm 
and the Perseus arm. It was first classified as a PN and is 
nowadays generally considered to be a compact H n region. 
To investigate the nature of Ml-78 further, we obtained long- 
slit, intermediate-resolution, optical spectroscopy with the ISIS 
spectrograph at the WHT. As a complement, we obtained long- 
slit, intermediate-resolution, NIR spectra using LIRIS, the NIR 
imager/spectrographer also installed at the WHT. 

Ml-78 is characterized by two main morphological zones, 
a bright arc to the SW and a blob of emission to the NE, which 
correspond to two peaks in the surface brightness distribution 
of the optical and NIR emission lines. In addition, each peak 
contains a featureless, very weak and unresolved stellar contin- 
uum that was detected in the optical. 

The electron densities calculated from the [S iij 6717/6731 
line ratio are found to be very high and similar in both zones, 
6700 and 6400 cm~^ for the arc and the blob. These are con- 
sistent with the densities determined from the NIR [Fe iii] line 
ratios. The electron temperature determined from the optical 
nebular [N ii] lines, on the other hand, is substantially different 
in the two zones, 10900 K in the SW arc and 13400 K in the NE 
blob. This last temperature is atypically high for a Galactic H ii 
region. The temperature determined from the NIR He i lines is 
somewhat lower, but this difference between Te obtained from 
coUisionally excited and recombination lines could be related 
to the presence of temperature fluctuations. We find as well lo- 
cal differences in the extinction. The spectrum of the NE blob 
shows a somewhat larger reddening coefficient than the SW 
arc, indicating a larger amount of internal dust in this zone that 
produces a higher Ay of about 9 mag. 

The most important result, however, is the confirmation of 
a nitrogen enrichment in Ml-78. This enrichment is stronger 
at the location of the NE blob. It seems that a process of local- 
ized chemical enrichment is taking place in the NE of Ml-78. 
This N enrichment in the blob is correlated with a defficiency 
in the O abundance and a - dubious - He enrichment. Such 
abundance pattern is typical of ejecta nebulae around evolved 
massive stars such as WR stars and LBV stars. The high exci- 
tation and electron temperature found in Ml-78 suggest that it 
is more likely to be a WR star: ring nebulae around WR stars 
show a wide range of nebular excitation and while LBV neb- 
ulae usually show lower than 10000 K and very faint or even 
unobservable [Om]. 

The spatial variations in the physical conditions and chem- 
ical abundances and the presence of more than one possible 
ionizing source indicates that Ml-78 is better described as a 
combination of a compact Hn region -I- ejecta. In fact, the 
Hei2112jum/Br7 line ratio, which is commonly used to esti- 
mate the spectral type of ionizing stars in H ii regions, indicates 
a hot {Teff > 40000 K) O star in flie SW arc. In confusion, Ml- 
78 seems to be ionized by a WR -i- hot O star system. Further 



observations able to observe directly the phosphere of these 
stars will be necessary to confirm this result. 

Finally, we detect H2 emission wich extends over a large 
(~ 30") area around the ionized nebula. The analysis of the 
NIR H2 lines, which includes estimates of the ortho-to-para ra- 
tios and excitation diagrams, indicates that the excitation mech- 
anism is UV fluorescence. We do not find evidence of shocks. 
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